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, Abstract. The string breaking model of particle production is extended in order to 

(N ■ 

help explain the transverse momentum distribution in elementary collisions. Inspired 
by an idea of Bialas', we treat the string using a collective coordinate approach. This 
leads to a chromo-electric field strength which fluctuates, and in turn implies that 
quarks are produced according to a thermal distribution. 
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1. Introduction 

O ' We study the production process of hadronic particles in elementary collisions based 
^ ' on the string breaking mechanism of QCD. The non-Abelian SU(3)-color gauge group 
^ ■ leads to confinement of the color electromagnetic field lines at scales of roughly 1 fm and 
above. Thus, a qq pair produced in, say, an e~^e~ collision will have a chromo-electric 
flux tube connecting the quark and antiquark. The tube increases in length as the 
Q ■ quarks separate, and eventually there is sufficient energy in the field to produce another 
qq pair. The rate of production is calculated via Schwinger's formula P, generalized 
to consider transverse particle momentum 0. This is often referred to as the string 
breaking mechanism of particle production. 
' Despite its successes, one troublesome detail of the QCD-string model has not 

been fully explained. The predicted transverse momentum distribution of quarks is 
Gaussian [2], while the observed spectrum of final state particles is thermal, or very 
nearly so Hj. Naive averaging over constituent quark transverse momenta cannot 
make an exponential distribution out of a Gaussian one. Further, it is not natural to 
suppose that there is sufficient rescattering for thermalization in systems as small as 
those in e+e~ and pp collisions. Hagedorn thus spoke before QCD was developed of 
'pre-established thermal equilibrium' While the string model itself has been refined 
several times, nothing like a thermal p± spectrum has been predicted by any of these 
refinements until recently. 

This discrepancy invites new ideas. Bialas showed that a random fluctuation of the 
string tension could produce particles in a thermal spectrum , if the appropriate initial 
distribution is used. He proposed that these fluctuations originate from the stochastic 
nature of the QCD vacuum. There has been further work on an extensions of this 
idea using dynamical fluctuations of the string in time [H]. This method seems like a 
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promising way to proceed, and the purpose of this paper is to improve the justification 
for these fluctuations, as well as to propose a possible origin. 

We separate the transverse and longitudinal dynamics of the string, introduce 
collective coordinates for the transverse dynamics of the string and quantize them. The 
reason that we will proceed this way is the following order of magnitude consideration: 
due to the very high momentum in the longitudinal direction (constituent quarks may 
have pl ~ lOOGe^ or more, in a 1 TeV pp collision) it is reasonable to treat this part of 
the dynamics classically. However, the typical transverse momentum of a produced 
qq pair will be on the order of the temperature (which we are seeking to justify) 
T ~ 160 MeV. The radius of the string is about 1 fm, approximately the same magnitude 
as the DeBroglie wavelength of the produced particles. 

Our approach is inspired by prior collective coordinate approaches in many-body 
quantum systems. One well-known example is the excitation spectrum of vibrations and 
rotations of large nuclei[7]. In the simplest liquid drop version of this model. Coulomb 
repulsion competes with nuclear surface tension. Quantization of the Hamiltonian 
resulting from these interactions leads to an excitation spectrum which can be verified 
experiment ally. 

In the string model, we once again have two competing energy considerations. The 
field lines will tend to spread in order to minimize their energy content. However, 
the vacuum confines color; this effect can be effectively reproduced by implementing 
in transverse direction the bag pressure and energy density dynamics. One might ask 
where the many-body quality arises in the string picture, since there is just a single 
qq pair anchoring the ends of the string. The answer is in the large number of virtual 
gluons which constitute the string QCD fields. We consider the string to be a collective 
excitation of this gluon sea. 

In section |2l we first look at the classical string picture, namely the energy balance 
and the Hamiltonian that results. Quantizing this Hamiltonian yields a harmonic 
oscillator-type wave equation for the transverse dimension of the string. We explore 
some of the properties of the solutions to this equation. Then, in section El we fold the 
resulting probability distribution for the string tension with the Schwinger formula for 
pair production [H |2] . This in fact generates an exponential p± spectrum. 

2. Wave Equation 

2.1. String Hamiltonian 

In preparation for quantization, we will first write down a Hamiltonian, focusing here 
on the energy density per unit length. We combine the contributions to the energy 
density from the chromo-electric field and the vacuum pressure to obtain the classical 
Hamiltonian density 

a = ^ElA + BA. (1) 
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The two terms are reminiscent of the usual kinetic and potential energies. Since the 
chromo-electric field is longitudinal, we will denote it by to avoid confusion with 
energy. 

Suppose the field lines are broken by a qq pair produced, each having color charge 
g. In the case where the field is longitudinal, we may simplify the non-Abelian Gauss' 
law down to its usual form: 

Jj2EAa-dA=^g. (2) 

The factor ^ arises comparing the elementary currents in QCD and QED: 

jH ^ I^T^'A^vP, ^ eijY^. (3) 

Note that we have taken El to be constant here. This is easily justified: under the 
constraint of a fixed fiux, the total energy is minimized by a uniform field. Due to the 
presence of the constraint equation Q, one ehminates the dependence on, say, E^, and 
obtains: 

.= |^ + SA (4) 

We can see the interplay between field and vacuum energy more clearly now. Minimizing 
a with respect to A gives 

yj^classical 9 /c\ 

~ 2V2B' ^ 

^classical ^ ^25, (6) 

and thus a classical string tension of: 

2.2. Position-momentum collective dynamics 

The string arises as a quasi particle state (presumably) in the full theory of QCD. Thus, 
its dynamics must be treated in an appropriate quantized way. However, as we noted 
earlier, the longitudinal evolution of the string should be safe to approximate classically, 
at least in our initial crude treatment. Furthermore, we assume the string is rigid. To 
be precise, there is no spatial bending of the string, nor spatial dependence in the field 
strength. 

There are two reasons that we make this assumption. First, though the energy may 
fiuctuate between bag and field energy, we suppose that the configuration of either is such 
that their separate energy contents are minimized. Second, treating the complexities of 
a "wavy" string is well beyond the scope of this paper. The chromo-electric field now 
contains only 1 degree of freedom, the longitudinal strength. Therefore, its conjugate 
coordinate, which will be the string cross section area or something related, will retain 
only 1 degree of freedom. Thus, it appears that the initially 2-dimensional transverse 
dynamics may in fact be successfully modeled by a 1-dimensional wave equation. Despite 
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the simplifications, it will still be a problem with a constraint, equation Q (Gauss' law), 
and an unusual looking Hamiltonian. 

To deal with the transverse dynamics as a quantum problem in an expedient 
manner, we make the following further simplification. As is evident the problem arising 
should still be closely related to the original one, though, admittedly, there is no absolute 
guarantee of complete isomorphic relation. We take the constraint equation (j^J and 
substitute it just once into the Hamiltonian, equation ((TJ, yielding: 

= ^gEL + BA. (8) 

We want to quantize this Hamiltonian in terms of Hermitian operators with the usual 
dimensions of distance and momentum. Since our degrees of freedom are A and E^, 
which dimensions are square of distance and momentum, we introduce 

x = VA (9) 

P=^L- (10) 

Here we have defined p to be the root of the classical energy density. Since the square 
root may take on positive or negative values, we will allow x, p e (— cxo, +oo). We make x 
andp into canonically conjugate operators, and so they satisfy the commutation relation 

[x,p\=ig (11) 

The justification for our choice of commutator will follow. Indeed, H„ now looks like 
the Hamiltonian of a harmonic oscillator: 

H, = ]^p' + Bx\ (12) 
Working in the p-representation leads to the wave equation: 

(^P' - Bg'^^i ^„ = (13) 
This is a very familiar equation, and its eigenvalues are: 

a„ = V^(r^ + ^). (14) 

Note that if the dimension of the problem is in fact not 1 but d, the | would have to be 
multiplied by c?, and there would be additional n's. 
As a matter of convenience, let us define 



(15) 



Y27r(2n + 1) 

Then we may express the eigensolutions of the problem as 

An 



Mp) = ^7^n(Vp74vrT2)e-^^/«-^o, (16) 
V -^o 

where 7i„ are Hermite polynomials, and An is a dimensionless normalization. 
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2. 3. Properties of the quantized string 

Let us now take stock of some of the features that have become evident in our model. 
First of all, we may calculate the expected cross sectional area and field strength; because 

= {Bx') = (17) 



we obtain 



{El) = — = V2Bi2n + l), 
9 



In both cases there is agreement between the ground state and the classical picture. 

Now, we want to see if the constraint, as well, is satisfied. It is reasonable to look 
at the product of the expectation values and see if Gauss' Law results. 

{EL){A) = ^-{2n + ir (20) 

By using this combination of operators we see that the ground state at least satisfies 
((21). This also clarifies our choice of commutator in equation (jllj) . Whether or not this 
excludes excited states of the string remains to be seen. For if we consider different 
combinations of expectation values that still appear to be related to the constraint, we 
attain different results. Consider, for example, 

(xp)(px) = |. (21) 

This should be related to the constraint, and it indicates that all excited string states 
may satisfy equation Finding the correct combination of operators to verify if 

the constraint equation Q is satisfied remains a difficult problem. However, these 
preliminary results strongly suggest that at least one, if not more, string states are 
admissible solutions. 

3. Thermal p± Spectrum 

We are now ready to investigate the spectrum of particles produced from this string. 
First, let us look at the probability distribution P{El). Taking the probability density 
if)^ and changing coordinates back to the field strength, we obtain 

Br, 



dP^E^) = -—^nl{^gEj8nTi)e-^^-/'-^odEL. (22) 

Bn is again a dimensionless normalization. Note that the factor {gEL)~^^^ arises from 
the change in measure. The singularity arising for even n is not troubling; it is an 

— 1/2 

artifact of the change in measure and never worse than E^ 

Next, we review the traditional particle production spectrum. The probability per 
unit 4- volume of a quark pair of momentum p± being produced is given by the Schwinger 
mechanism, as expounded by Casher, Neuberger and Nussinov[2|: 

dP{p^) = ^e-2-</^^MV- (23) 
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This is in fact only the first term in an infinite sum, but the approximation is appropriate 
for our current level of precision. We consider each string breaking to be, essentially, a 
measurement of El. For now, let us suppose that the strings are formed in the ground 
state, with string tension (Tq. Thus, over several string breakings, the observed spectrum 
of produced quarks is obtained by folding equations and over El- 

Performing this folding results in the quark transverse production probability per 
unit 4- volume, 

dP^u.AP±) = ^(1 + m^/To)e-"^-/^«rfV. (24) 
9 

This agrees with previous results up to a linear prefactor[5J |^. As a final step, consider 
the composite hadron spectrum. We will obtain this by first rewriting the distribution 
in terms of the transverse mass (fortunately, it is already more or less in such a form), 
and then folding two distributions together, corresponding to quark-antiquark or quark- 
diquark pairings. The final hadron spectrum thus obtained is 

C^Phadrons(m^) ^ / _ 6 ^ ^ ^ ^^^^ 



dm± \ 5 



-Xl (3X2 + 2X3) +Xl + Xi + ) e-^^ (25) 



where 

X^^ ( , , , , , _ 

nJ \ToJ ' ^ To 

and rui and m2 are the constituent quark/diquark/antiquark masses. Note that the 
largest root of the prefactor is m_|_ = mi + m2, and thus the probability is always 
positive, as expected. This can be seen trivially in the case uii = 1712 = 0. The specific 
relationship ()15p is the same as in the aforementioned works, though the pre-factor 
shifts slightly the fitted temperature. Making use of a fairly standard value for the 
string constant, ctq ~ 0.9 GeV/fm, gives Tq ~ 165 MeV. 



4. Conclusions 



We have considered the string between a qq pair as a quasiparticle excitation of the gluon 
field. Thus, we justified introducing collective quantization of transverse dynamics. 
This implies that the internal chromo-electric field is not constant, but rather when 
the string breaks and the field is observed, it takes on values according to a quantum 
probability distribution. Folding this distribution with the traditional Schwinger- 
Casher-Neuberger-Nussinov result leads to an intrinsically thermal p± spectrum. The 
emerging phenomenology is quite satisfactory; e.g. it produces a thermal spectrum. 
However, a more complete theory will certainly lead to further insights including 
refinement of the temperature parameter. 
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